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Abstract 

Two implicit periodic structures in the solution of sinh-Gordon thermodynamic Bethe 
ansatz equation are considered. The analytic structure of the solution as a function of 
complex 9 is studied to some extent both analytically and numerically. The results make 
a hint how the CFT integrable structures can be relevant in the sinh-Gordon and staircase 
models. More motivations are figured out for subsequent studies of the massless sinh-Gordon 
(i.e. Liouville) TBA equation. 

1. Sinh-Gordon model 

It seems that from the recent developments of the string theory there are some persistent 
requests about better understanding of the two-dimensional sigma-models with non-compact 
(in particular singular) target spaces. Physical properties of such models are expected to be 
quite different from these of the better studied compact sigma-models. It is therefore a chal- 
lenge for the two-dimensional integrable field theory community to reveal the corresponding 
peculiarities and new features. 

The two-dimensional sinh-Gordon model (ShG) is defined by the (Euclidean) action 
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ShG 



— (d a 0) 2 + 2/icosh(260) 
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d x 



I believe that this model can be considered to say as "a model" of (suitably perturbed) 
non-compact sigma-models, in the sense that its properties are considerably different from 
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those of the perturbed rational conformal field theories, these differences being sometimes 
quite similar to these between non-compact and compact sigma-models. 

In (|1.1|) /i is a dimensional (of dimension fi ~ [mass] 2+2fe2 ) coupling constant which deter- 
mines the scale in the model and 6 is the dimensionless ShG parameter. For the beginning 
I suppose it to be real and non-negative. The case 6 = turns to be somewhat singular for 
the subsequent considerations and here the study is restricted to positive values < 6 < oo 
only. It is also convenient to use different parameters 

Q = b+l/b (1.2) 

and 

instead of b. 

The perturbing operators exp(±26</>) in ( |1.1|) have negative dimension A = —b 2 . To make 
the coupling fi a strict sense we need to fix a normalization of these operators. Here they 
are implied to be normal ordered w.r.t. the massless (unperturbed) vacuum in the way that 
in unperturbed theory 
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The model is known to be integrable and can be solved for many important characteris- 
tics. In particular its factorized scattering theory is known since long The spectrum 
consists of only one neutral particle B(9) subject to a factorized scattering with two-particle 
amplitude 

smh V + % sm 7rp 

With the normalization (p..4|) the mass m of this particle is related to the scale parameter fi 
as § 

^ J {b \ = [mZ(p)] 2+2b2 (1.6) 



r(i-6 2 ) 



where 



z (p ) = ^(i- P )^r(l^)r(|) (1.7) 

Notice that the scattering theory is invariant under the (week-strong coupling) duality 
transformation b — > 1/6 which brings p — ► 1 — p. This means that the physical content of the 
model remains unchanged up to the overall mass scale. Since the combination ( |1.7| ) is again 
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invariant under p — > 1 — p, the mass scale also remains unchanged if the coupling constant 
/i is simultaneously substituted by the "dual" coupling constant fi related to fi as follows 

T(b 2 ) \ 1/b ( Y(\lb 2 ) \ b 

Therefore the sinh-Gordon model is completely invariant under duality b — ► 1/6; // — > jl. 
Due to this symmetry it is sufficient to consider only the region 0<6 2 <lor0<p<l/2. 

The infinite volume bulk vacuum energy of the model is also known exactly [|J. In terms 
of the particle mass m it is given by the following apparently self-dual expression 

8 sinn up 

2. TBA equation 

Contrary to the on-mass-shell data of ShG quoted above, the off-mass-shell characteristics 
such as the correlation functions (with the exception of the vacuum expectation values of 
some local fields 0,|6| and their matrix elements between the asymptotic states J7j) are not 
known exactly 

Some progress can be made for the finite size effects where the problem is reduced to a 
non-linear integral equation known as the thermodynamic Bethe ansatz (TBA) one. Namely, 
consider the ground state energy E(R) of the finite size ShG model placed on a circle of finite 
circumference R. In the TBA framework it appears as 

E(R) = -— [ cosh 6 log (1 + e- £{e) ) d6 (2.1) 
2n J 

It is convenient to introduce also the independent effective central charge c e g(i?) as 

c eS (R) = --E(R) (2.2) 

7T 



In (2.1) e(6) is the solution of the TBA equation 

mR cosh 9 = e + <p * log (l + e~ £{9) ) (2.3) 
where * denotes the ^-convolution. The kernel if (8) is related to the ShG scattering data 

(PD 

<p(e) = -~*>gS(e) = - 4 f7 coshg (2 .4) 
2% dO 2n cosh 29 — cos 2np 

The Fourier transform of the kernel reads 



cosh • 



ip(u) = I e»»<p{9)d0 = 4jt (2.5) 

cosh — 
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where the parameter 

a = 1 - 2p (2.6) 

is simply reflected a — > —a under the dualily p — > 1 — p and therefore can be taken non- 
negative < a < 1. 

The following conclusion are readily made from the structure of the integral equation 

1. Function is even = e(—9) and analytic in the strip |Im#| < tt/2 — na/2. At 
Ke9 —>■ oo in this strip it has the asymptotic e(9) ~ mRe e /2. Therefore the function 

= exp(-e(0)) (2.7) 

is analyic and non-zero in this strip and at Re 9 — > oo behaves as 

~ exp(-mRe 6 /2) (2.8) 

The asymptotic at Re# — > — oo is related to (|2.8| ) by the symmetry Y(8) = Y{—9). Let us 
define another even function 



X{6) = exp 



mR , n /" log(l + Y(B')) dff 
■cosh 6>+ ' 



2sin7rp J cosh(6 ) — 6') 2n 



(2.9) 



which is obviously analytic and non-zero in the strip | Im 6 \ < it/ 2 and at Re 6 — > oo in this 
strip 

From the TBA equation it follows that 

X(6 + ian/2)X(6-iair/2) = Y(6) (2.11) 

2. In fact on the real axis Y{9) is real and positive and therefore we expect a strip 
| Im0| < e with some finite e > where 1 + Y(9) ^ 0. Therefore the analyticity condition 
for X{9) can be extended to the strip | \m.6\ < n/2 + e. This is enough to prove the relation 

X(6 + m/2)X(6-iir/2) = l + Y(9) (2.12) 

The functional equation 

X (9 + ztt/2) X (9 - in/2) = 1 + X (6 + ian/2) X (9 - ian/2) (2.13) 

follows. This relation allows to extend the original analyticity strip | Im#| <7r/2 + eto the 
strip | Im#| < 3n/2 and, as we'll see before long, to the whole comlex plane of 9 so that 
X{9) is an entire function of 9. Notice that from (|2.13 ) it follows that the asymptotic ( |2.10|) 
holds in the larger strip | Im#| < n. The asymptotic outside this strip is more complicated. 
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3. As a consequence, Y(9) is also entire function of 9 and satisfyes the following functional 
relation 

Y (9 + in/2) Y(9- in/2) = (1 + Y (9 + ian/2)) (1 + Y (9 - ian/2)) (2.14) 

The last equation is very similar to the functional relations appearing in the TBA study 
of the integrable perturbed rational conformal field theories (the so called F-systems) . Typi- 
cally such F-systems imply a periodicity of the F-functions in 9 with some imaginary period 
related to the scale dimension A of the perturbing operator (see e.g., ||[9||). This periodicity 
in order entails special "perturbative" structure of the short distance R — > behavior of the 
ground state energy E(R). Namely, up to one exceptional term, it is a regular expansion in 
powers of R 2 ~ 2A 

00 

E{R) = -S vac R --^Y, c nR (2 ~ 2A)n ! R - (2.15) 

n=0 

where £ vac is the infinite volume vacuum energy of the model. Unlike this typical situation, 
the F-system ( j2.14|) or the X-system ( [2.131 ) does not imply any apparent periodic structure 
of Y{9) in 9. As a manifestation of this peculiarity, the R — > behavior of E(R) is different 
from (|2.15|) and includes softer logarithmic corrections 0] 

E(R) = -SR-^(l- ^ g2 ~ R P) + O (log" 3 it!)) ; R - (2.16) 

The purpose of the next two sections is to reveal two hidden periodic structures (with different 
periods) of the F-system ( p,14|) . 



3. Discrete Liouville equation 

In this section I discuss the following two-dimensional non-linear finite-difference equation 
for the function X(u, v) 



X(u + 1, v)X(u -l,v) = l + X(u, v + l)X(u, v-1) 



(3.1) 



which is apparently resemblent of the functional A-system ( p,13|) . In the next section we'll 
see how some of the results for ( |3.1|) can be specialized to our TBA problem. Equation 
( |3.1D is a particular case of Hirota difference equation ||11|| . The constructions of this section 
can be found in |12| (see also [J3]) where more general difference system is analysed. They 



appeared also in a quite close context in [11 . 

Equation (|3.1|) can be considered as a discretisation of the hyperbolic Liouville equation 



ftp - d lv = ~ e 



2y 



(3.2) 
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Indeed, let X(u,v) = exp(—<p(u,v)) and let ip(u,v) be large and negative. Then eq.(|3.1|) is 
approximated by 

ip{u + 1, v) + ip(u — l,v) — <f(u, v + 1) — ip(u, v — 1) = — exp(<p(u, v + 1) + f(u, v — 1)) 

(3.3) 

In the long- wave limit this is reduced to ( p.2|) . As we'll see below, eq. (|3.1|) is in many respects 
very similar to the Liouville equation. It seems quite natural to call it the discrete Liouville 



equation 15 



Let me remind well known construction of a local solution to eq.( |3.2| ). 
la. It is convenient to use the light cone variables x + = u + v and x~ = u — v , so that 
d + = (d u + d v )/2; 9_ = (d u - d v )/2 and fO) reads 



Ad+d^ = -e 2ip (3.4) 



Let (p be a local solution of (|3.4j) . Define 

t = -{d +V f + d 2 + cp (3.5) 

t = -(d^f-d 2 ^ 

As a consequence of eq.(p.4|) we have 

8_t = d + t = (3.6) 

so that t(x + ) and t(x~) are respectively functions of only x + and x~. 
2a. Field X = exp (—(f ) satisfies two linear differential equations 

(d 2 + + t(x + ))X(u,v) = ; (d 2 _+t(x-))X( Ul v) =0 (3.7) 

3a. Let Q±(x) and Q±(x) be linearly independent solutions to the ordinary differential 
equations 

(dl + t(x))Q±(x) = (3.8) 
(d 2 x + t(x))Q±(x) = 

normalized in the way that 

d x Q + (x)Q-(x) - Q + (x)d x Q4x) = 1/2 (3.9) 
d x Q+(x)Q4x) - Q + (x)d x Q4x) = 1/2 

A local solution of (|3.4| ) can be constructed as 

exp(-<p{u,v)) = Q + (x + )Q + (x-) + Q^{x + )Q^{x-) (3.10) 
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4a. Introduce the functions 



a(x) -W) (3 ' n) 

The solution (|3.10|) can be rewritten as 

exp¥v)= (^F« (3 - 12) 

Notice also that in terms of F and G 

t(x) = -2{F(x),x} ; = -2{G(x),x} (3.13) 

where 

= y " 2 (7) ( 3 - 14 ) 

is the Schwarz derivative. 

Now let us turn to the discrete Liouville equation (|3.1|) . Defining, similarly to eq.( |2.11|) , 

Y(u,v) = X(u,v + l)X(u,v- 1) (3.15) 
1 + F (w, v) = X(u + 1, v)X(u - 1, u) 

we arrive at the finite difference equation analogous to ( |2.14|) 

Y{u + 1, v)Y(u - 1, v) - (1 + F(m, u + 1))(1 + y(u, u - 1)) (3.16) 

lb. Introduce 

Tf v X(^ + l,t; + l)+X( W -l,t;-l) 

v) = ( 7) 

X(u, v) 

As a consequence of eq.( |3.1D we obtain 

T(u + l,v-l)=T(u,v) (3.18) 
f(u + l,v + l) =f(u,v) 

so that T = T(u + f) and T = T(w — u). In the continuous limit these objects are related to 
(p|) as T(u) = 2 - 4t(u) + ...; f{u) = 2- At(u) + ... 0. 
2b. Eqs.( |3.l7D can be rewritten as (similarly to ( |3.7| )) 

X(u + 1, v + 1) + X{u - 1, v - 1) = T(u + v)X{u, v) (3.19) 
X(u + l,v-l) + X(u -l,v + l)=f(u- v)X(u, v) 
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3b. Let Q±(u) and Q±(u) be linearly independent solutions of the second order finite 
difference equations 

Q ± (u + 2) + Q±(u-2) =T(u)Q±(u) (3.20) 
Q ± (u + 2) + Q ± (u - 2) = f(u)Q ± (u) 

normalized by the "quantum Wronskians" 

Q + (u + l)Q-(u - 1) - Q + (u - + f ) = f (3.21) 
Q + (u + l)Q-{u - 1) - Q + (u - l)Q-{u + 1) = 1 

Then it is verified that 

X(u, v) = Q + {u + v)Q+(u -v) + Q-(u + v)Q-(u - v) (3.22) 



is a local solution of the discrete Liouville equaiton (|3[ 
4b. Introduce the functions 

F(u) = ^l; GH = -lM (3.23) 

QA U ) Q+(u) 

which can be used to present the local solution ( |3.22| ) in the form 

(F(u + v + l)-G(u-v-l))(F(u + v-l)-G(u-v + l)) 
1 ' ' (F(u + v + l)-F(u + v-l))(G(u-v + l)-G(u-v-l)) { ' 
_ (F(u + v + l)-G(u-v + l))(F(u + v -l)-G(u-v- 1)) 
+ (u, v) - ^ + V + V) _ F ( u + V _ i))(G(u -v + l)-G(u-v-l)) 

Let me mention also the discrete analogue of the Schwarzian derivative ( |3.13| ) 

T(u + l)T ( u - 1) = \ F <" + 3 > ' F '» - W + 1) - f (« - «» ,3.25) 
v 7 v ; (F(u + 3) -F(u+ l))(F(u- 1) -F(u-3)) v ; 

4. Application to TBA 

The above constructions for the discrete Liouville equation can be translated for the ShG 
X-system (|2.13|) we are interested in. Let us require the following periodicity condition for 
X(u,v) in (Q 

X(u + a,v + l) =X{u,v) (4.1) 



with some parameter a (at this point we start to diverge from the lines of [Q). With this 
periodicity eq.( [3.1|) reads 

X(u + 1, v)X(u - 1, v) = 1 + X(u + a, v)X(u - a, v) (4.2) 
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Here v can be considered as a parameter. Suppressing this redundant dependence and 
rescaling u as 



9 = inu/2 (4.3) 

we are back to the ShG X system in the form ( [2.13|) . 

From X{9) the two functions T(9) and T(9) are readily restored as 

X(fl + m(l-a)/2)+X(g-m(l-q)/2) (A ^ 

= ^ (4.4) 

- X (g + i7r(l + a)/2) + X (g - »tt(1 + q)/2) 

TW = W) 

The "holomorphic" property (|3.18| ) is translated to the following periodicity of these func- 
tions 

T{9 + m{l + a)/2) = T{9) (4.5) 
f(6 + t7i(l - a)/ 2) = f{9) 

Notice that the period Z7r/(1 + b 2 ) of T corresponds to the negative dimension A = —b 2 
of the perturbing operator in ( |1 . 1|) . As it can be anticipated from the self-duality of ShG, 
the second period mb 2 j(\ + b 2 ) of T is related to the dimension A = — b~ 2 of the "dual" 
exponentials exp(±20/6). 

As it is discussed in sect. 2 X(9) is analytic and non-zero in the strip | lm9\ < it/2 and 
analytic in the larger strip | Im#| < 37r/2. Therefore T(6) and T(6) are analytic in the strip 
|Im#| < 7r/2 and by periodicity ( [4.5[ ) are entire functions of 9. It follows from ( |4.4| ) that 
X(8) is an entire function of 9 too. 

The asymptotics at Re 9 — > oo follow from (|2.10 ) and (|4.4| ) 

(4.6) 



TV) ~ exp ( m ^ eXP 4 ( l; ( ^ ( / 2) ~ P)/2) ) - the strip < Imtf < .(1 + «)/2 
f{6) ~ exp f ^gj^jWg)^ in the strip o < ke < 7r(1 _ a)/2 



4sin(7rp/2) / 
The real axis Im # = is a Stokes line and here 



T{9) ~ 2 exp ( ^ tan(7rp/2)e^ cos exp 5 ) (4.7) 

, mi? . . , n\ ( mR 
T(9) ~ 2 exp — — cot(7rp/2)e e cos — — exp# 



The lie 9 — > oo asymptotic in the whole plane of 9 is restored from the periodicity 
Following (|4.7| ) both T and T have infinite number of zeroes on the real axis located at 
9 = ±6 n , n = 1, 2, . . . , oo with 9 n ~ log(27rn/mi?) + 0(l/n) at n — > oo. The half-period 
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shifted functions T{6 + m{\ — p)/2) and T(9 + iirp/2) are also real at real 9 and at 9 — > oo 
behave as 

r(e + i .( 1 - p)/2 ).exp( i ^ M e») (4.8) 



T{9 + inp/2) ~ exp 



mi? 



4sin(7rp/2) 

Analytic properties of T and T allow the following convergent expansions 

oo 

T{6) = J2 Tn exp(2nQ&0) (4.9) 

— oo 
oo 

f{0) = Y,T n exp(2nQ9/b) 

— oo 

with real T n and T n . From the symmetry of the original TBA equations T and T both even 
functions of 9 so that in our present case T n = T_ n and T n = T„„. One can read-off the 
leading large n behavior of the coefficients T n and T n from the asymptotics ( ^4.6| ) 

Tn „ (_)» t M n -*-V3 L mR \ 2QbU ( 4.io) 



7r \ 8irbQ cos(irp/2) 

T n ~ (-)»J9- n -*»W* (e hmR ..Y Qn/b 




'8ttQ sin(7rp/2) J 

So far the constructions were explicitely based on the integral equation The rest 

of the section is more speculative. Suppose that for T{9) and T{6) constructed as in (|4.4| ) 
we can find Q±{9) and Q±{9) which solve 

Q ± (9 + in) + Q ± {9 - m) = T{9)Q ± {9) (4.11) 
&t (0 + in) + Q±(0 - itt) = T{9)Q ± {9) 

and are the "Bloch waves" with respect to the periods of T and T respectively 

Q± (9 + itt(1 + o)/2) = exp (±2ivrP/Q) Q±(0) (4.12) 
Q± (0 + in(l - o)/2) = exp (±2inP/Q) Q±(6) 

with some Floquet index P. Let them be normalized by the quantum Wronskians 

Q + (9 + in/2)Q_(9 - in j2) - Q + (9 - nr/2)Q_(0 + ztt/2) = 1 (4.13) 
Q + (9 + z7r/ 2 )<2-(# - ztt/2) - Q + (9 - z7r/ 2 )<2-(# + W 2 ) = 1 

Then formally 

X(9) = Q + (9)Q + (9) + Q-(0)Q-(0) (4.14) 
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solves eqs.(PO) as well as the X-system ( |2.13| ). 



Unfortunately at present I know no effective means to construct these Q-functions. More- 
over, there are serious doubts that the objects satisfying both ( [4.11 ) and ( |4.12|) can be con- 



structed in any sense, at least at rational values of b 2 . I hope to say something more definite 
on this point in close future. 

5. Large Re# asymptotics 

Let me comment a little more about the Re# — > oo asymptotics (with lm.6 fixed) of the 
function X{9) in the whole complex plane. In principle it can be restored from the asymptotic 
( [2.1(J| ) in the strip | Im#| < 7r/2 using the functional relation ( |2.13| ) or, more conveniently, 
the relations (|4.4j) together with the asymptotics ( |4.6|) . The asymptotics is always of the 
form 

X{6) ~ exp(A(Im#)exp(Refl)) (5.1) 

with some complex function A(rj) of real variable rj = lm.6. Apparently Re A(rj) controls the 
rate of growth of the absolute value of X. At |?7| < tt /2 we have 

MV) = --^—^ (5.2) 



It follows from (|4.4j ) that A{r\) satisfies two functional relations 

A(r, + ,(1 - a) ,2)) = (- 4M + I Re( f "> > + 1(1 " ° )/2) (5.3) 

w v " " 1 A(t7 + 7t(1 -a) /2) otherwise v ; 



A(r) + ir(l + a)/2)) 



fA(r]) + a(rj) if Re(A(rj) + a(r])) > ReA(r/ + tt(1 + a)/2) 
\ A(r] + 7r(l + a)/2) otherwise 



where the functions a(j]) and 17(77) control the asymptotics of T(6) and T{9). They are 
defined as 

<?{V) = A mR /0 , exp(^ - £tt(1 - p)/2) at < 77 < (1 + a)/2 (5.4) 
4cos(7rp/2) 

^ = /1 • /o\ ex P( ir 7 - inp/2) at < 77 < (1 - a)/2 

4 sm(7rp/2j 

and continued outside these regions periodically as a(rj + 7r(l — a)/2) = o"(r/) and a(r] + 
7r(l + a)/2) = (j{rf). Both cr(r/) and &{rj) jump by —imR/2 at 77 = tt(1 + a)n/2, n e Z and 
77 = 7r(l + a)n/2, n G Z respectively. This corresponds to the limiting density of zeroes 
prescribed by (fOj). 

A common solution to (|5.3| ) exists. Generally the solution is discontinuous at all values 
77 = ±7r(m(l + a)/2 + n(l — a)/2) with arbitrary positive integers n and m. At each 
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Figure 1: Real and imaginary parts of A{rj)/mR at b 2 



such point the imaginary part ImA jumps down by —mR/2 indicating an asymptotic line 
of accumulation of zeroes of X(6), the asymptotic density being the same as that of the 
functions T and T ( |4.7| ). The real part of A at these points is continuous itself but has 
discontinuities in the first derivative. The first Stokes line appears at 77 = 71. 

At large 77 the structure is qualitetively different dependent on the arithmetic nature of 
parameter b 2 . If it is a rational number the periods of T and T are commensurable. Some 
of the discontinuities merge forming multiple jumps in the imaginary part. The solution 
A(6) bears a regular "quasiperiodic" structure with the common period of T and T . For 
irrational b 2 the periods are incommensurable and as 77 —>■ 00 the singularities are more and 
more dense, the solution having quite irregular behavior. 

In fig.l real and imaginary parts of A (77) are plotted for the simplest (self-dual) case 
6=1. Both periods are equal to 7r/2. The discontinuities are located at 77 = ±7r(n + l)/2, 
n = 1, 2, . . . where the imaginary part jumps by —mRn/2. Contrary to the location of zeroes 
of T, zeroes of X cannot lie exactly on the lines lm6 =const, at least for the first line 77 = ir. 
Indeed, in the case b 2 = 1 it follows from the functional relation that at real 9 

\X(6 + m)\ 2 = X 2 (6)+T 2 (6) (5.5) 

which is strictly positive. 

Another rational situation b 2 = 1/2 corresponds to the periods 7r/3 and 27r/3. Function 
A{rj) is plotted in fig. 2. The structure is again quite regular, the discontinuities occuring at 
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Figure 2: Real and imaginary parts of A(rj)/mR at b 2 



1/2. 



7] = tt + mr/3, n = 1, 2, . . . , the first two discontinuities in ImA being —mR/2, next two are 
twice of this amount, then next two trice, etc. 

With irrational b 2 the picture is far less regular. To illustrate what happens when b 2 
slightly deviates from a simple rational number, in fig.3 we plot A{j]) for b 2 = 0.8086 . . . 
which is reasonably close to 1. Comparing with fig.l we see that the first discontinuity at 
7] = ir remains basically the same while the second (double) discontinuity at rj = 3tc/2 splits 
in two simple ones, the third (triple) splits in three simple jumps etc. At some point these 
splitted groups come to overlap and the picture turns irregular. 

6. Staircase situation 

The staircase model 0] is a formal analytic continuation of ShG to complex values of the 
parameter b such that b~ l = b*. Although the physical content of this continuation is 
not completely clear from the field theory point of view, the TBA equation (|2.3|) remains 
completely sensible and this continuation of ( |2.3j ) can be studied on its own footing. The 
effective central charge fl2.2|) is still real and develops sometimes quite intriguing patterns 
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Figure 3: Real and imaginary parts of A(rj)/mR at b 2 = 0.8086 



with real < 6q < oo. Parameter p of (|1.3 ) now reads 



1 , i6 
p= 9 + ~T 

Z 7T 



while 1 — p = p* and a defined in Q2.6|) is purely imaginary a 
is real and reads 



<p{0) 



1 



1 



1 



2tt Vcosh(^ + 0„) cosh(# - # 



with the Fourier transform 



cos(uj9o) 
cosh(7ra;/2) 



(6.2) 

-2i9 /ir. The TBA kernel 
(6.3) 

(6.4) 



After these substitution the integral equation ( |2.3| ) determines real-analytic functions 
e(6), Y{9) and, through (|2.9|) , a real- analytic and symmetric with the asymptotic 

behavior at Re^ oo in the strip — 7r/2 < \m6 < 7r/2 



X(6>) ~ exp 



mR 



4 cosh #o 



exp 9 



(6.5) 
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The functional equation ( |2.13| ) reads now 

X (6 + ztt/2) X (6 - m/2) = I + X (6 + 6 ) X (6 - O ) (6.6) 

All the considerations of sectJJ can be repeated literally. Function X{9) is still entire as well 
as T(Q) and T{6) defined in eg. ([14]) and the asymptotic (|6.5|) can be extended to the strip 
—7i < lm8 < tt. The difference is that the periods of T (9) and T{9) 

T(6 + T ) =T(0) ; r = i7r(l + a)/2 = i7r/2 + o (6.7) 
f (9 + r) = f (9) ; f = -a) j2 = in/2 - 6 

are now complex r = — f*. Functions T(6) = T (—9) and T{9) = T(—6) are still symmetric 
but no more real analytic. Instead 



Expansions similar to (|4.9|) 



T*{9) = f{9*) (6.8) 

oo 

T{9) = T n exp(2iim9/r) (6.9) 

— oo 
oo 

f{9) = y *Tf n exp(2iim9/f) 

— oo 

are convergent and T n = T_ n ; T n = T_ n . Instead of being real as in the real b case, these 
coefficients are complex conjug cite J. ^ — J. . 

In the strip < Im6 < tt/2 the following Re# — > oo asymptotics holds for T and T 

T(S) ~ jr(g + f)/X( 9 ) ~ exp ( '" fi(1 4 ^2 ( ^ 0)) ^) (6.10) 
f (9) ~ X(9 + r)/X(g) ~ exp ( -= 



4 cosh 

The asymptotic changes at the Stokes line along the real axis where 

rp/n\ o ( mR ^ + * Sillh A f mR \ , r i i 

T(#) ~ 2 exp — e cos e (6.11 j 



4 cosh 6»o / V 4 

rp(n\ o /mi?(l -isinh^o) A /mi? 

T(6>) ~ 2 exp — e cos 

w V. 4 cosh 6» / 

and we observe again an infinite sequence of zeroes accumulating at infinity, the density being 
the same as in the real b case of sect.[5| In the strips irn/2 <Im# < 7r(n+l)/2 (with arbitrary 



15 



0.1 



-0.1 



Reo 




i i i i 


/ Ima 


r\/n 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

(a) 




Figure 4: Real and imaginary parts of a(r])/mR (a) and a(r))/mR (b). The staircase example 
for #o = 1. 



integer n) the Re# — ► oo asymptotics follow from the periodicity 
the lines Im# = inn/ 2 (with any integer n) we'll have 



In particular, along 



T{9) 
f(6) ~ 2 exp 



0(>X]) ( m - R ( 1 +» sinhg o) c Ree-ng 



4 cosh 6> 
mR(l — i sinh 0, 



cos 



4 cosh 6*n 



mR 
< 

4 

mi? 



,Re9-n6 



(6.12) 



,Re e+n8 



Functions cr(^) and «r(r/) which control the asymptotics of T and T at Re# — ► oo 

T{9) ~exp (a(Im^)e Re0 ) 
f(0) ~exp (&{Im6)e Red ) 



(6.13) 



are plotted in fig.[| for the case Oq = 1. It's enough to present them for n > since 
cr(-rj) = &*{rj). The imaginary part of (r{rj) jumps at the points rj = nn /2 by the amount 
—mRexp(—n8 )/2, in accord with the density of zeroes predicted by ( |S.12j ). Unlike the 
previously considered case of real b, in the staircase situation the zeroes of T and T are 
not located exactly at the lines lm8 = inn/ 2 but slightly shifted in the imaginary direc- 
tion (we'll observe this deviation numerically in the next section) and approach these lines 
asymptoticlally as Re 6* — > oo. 
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Figure 5: A staircase example of A(r])/mR at 9q = 1. 



Notice also that e.g. T(8) is a single valued function of £ = exp(2z7r^/r). In the complex 
plane of this variable the asymptotic lines of accumulation of zeroes lm.8 = 0; Re# — > ±oo 
are parts of the spiral |£| = exp(7r arg£/2# ) near which zeroes become dense at |£| — > oo or 
|£| — ► 0, the density growing as |£| ± ( 1 / 4+0 o/ 7r ) respectively. Therefore the large (or small) |£| 
asymptotics of T(£) at fixed arg£ is rather complicated. 

Large Re# asymtotic of X(8) at fixed rj = lm6 is controlled by the function A(rf) (see 
eq. (|5.1|) ). An example corresponding to the case #0 = 1 is presented in fig.|[ At the 
points f] = ±7r(n + l)/2, n = 1,2,... imaginary part of A(rj) has discontinuities equal 
to — m.Rsinh(n#o)/(2sinh#o)- These amounts determine the asymptotic density of zeroes of 
X(9) along these lines. 

7. Numerics 



Integral equation ( |2.3| ) can be easily solved numerically e.g., by iterations. The iterations 
happen to be well convergent (the convergence is somewhat slower if R approaches to 
or the parameter p is taken very small). In the strip | lm$\ < 7r/2 function X{6) can be 
computed using the integral representation ( |2.9|) . This allows to continue X{6) to the whole 
complex plane iterating the relation ( 2.13|) (in fact at large Im^ it is more convenient to 
evaluate first T{6) inside its period and then use flO]) ). In the rest of this section we will 
use the logarithmic scale parameter x = \og(mR/2) instead of R. 
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Figure 6: Typical bell-shaped patterns of X(9) on the real axis (b 2 = 1). 



1. Self-dual point b 2 = 1. 

In fig.|| several examples of function X(9) on the real axis of 9 are plotted for different 
values of x. Function is typically bell-shaped. As x becomes large negative the width of the 
support of the bell as well as its height grow proportionally to —x. No plateau typical for 
perturbed rational CFT's is developed. 

Few samples of T(9) (which is the same as T{9) at the self-dual point) are presented 
in fig.|7|. At x negative and large enough, T develops a plateau in the "central region" 
x < 9 < —x of the height which approaches slowly to 2 as — x grows. We'll comment more 
about this approach below. Outside the central region it starts to oscillate with growing 
amplitude and friquency. Approach to the asymptotic (|4.7|) is very fast. 

Due the the symmetry T{9) = T{—9) this function is real on the imaginary axis too. A 
couple of examples are plotted in fig.||. At x essentially negative, when the plateau is well 
developped in the central region, the mean value T (see eq.( |4.9|) ) is very close to the plateau 
height, the oscillations around (determined mainly by 7\) being very small (Ti ~ R 4 , see 
eq.flpD below). 

Function T{9) is real also at the half-period line lm.9 = 7r/4. Again there is a plateau 
region (at large negative x) of the same height as on the real axis. Then T{9 + z7r/4) remains 
positive and grows following the asymptotic (|4.8|) (see fig|§). Numerical comuputations in 
the whole period strip ^Im^ < 7r/2 show no sign of other zeros then those on the real axis. 

In fig.[Ty few examples of |A(6*)| in the complex plane at different values of Re 9 are 
plotted vs. Im^ (for x = 0). The specific values Re9 = 1.22 and Re9 = 2.25 are chosen 
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Figure 7: Function T(6) on the real axis (b 2 = 1). 
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Figure 8: Function T(9) along the imaginary axis (b 2 = 1). 
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Figure 11: First 10 zeros of X{9) located near the line lm.6 = tt. An example for x = and 
b 2 = 1. 



close to the positions of the first two zeros of T(9) on the real axis. The deeps near lm.6 = tt 
and lm.8 = 3n/2 indicate a presence of zeros of X(9) nearby. 

More precise positions of zeros of X{9) near the line Im^ = 7r (for the same value x — 0) 
are examplified in fig.[ll]. In fact all these zeros are inside the strip | \m6\ < tt. Only the 
first zero deviates noticably from the line ImO = tt. The imaginary parts of next zeros are 
already very close to tt and tend to this value very fast. 

2. Rational points. As an example of a rational point we take the simplest case 
b 2 = 1/2. The periods of T{6) and T(9) are commensurable and equal to 2m /3 and 27r/3 
respectively. In fact in this case there is no need to study separately T and T since, as it is 
readily derived from ( p. 13 ), they are bound up by the relation 



T(0) = T(6)T(6 + in/3) - 2 



(7.1) 



It should be noted that similar finite degree functional relations between T and T exist 
for any rational b 2 . For example, at b 2 = 1/3 the periods of T and T are 3iir/4 and in/A 
respectively and 



f{6) = T(6)T(0 + m/A)T(d + m/2) - T{6) - T{6 + zvr/4) - T(0 + in/2) 



(7.2) 



Numerical patterns of T (6), T{9) and X(9) are essentially the same as for b 2 = 1. I'd only 
show few plots of |-X"(0)| in the complex plane along the lines lm.6 = tt, Ait/3, 5tt/3, 2tt, Itt/3 
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etc. to illustrate the mechanism of multiplication of the zeros density in the asymptotics 
Re 9 — ► oo as required by the prediction of fig.|2|. In fig.[12| is plotted at Re9 in the 

vicinity of the real position of the first zero in T(9) at 9 = 1.2241 . . . (the case x = is taken 
as an example). At lm.9 = n and 47r/3 the picture indicates simple zeros located closely to 
this point in Re# and slightly displaced in the imaginary direction. At Im^ = 5tt/3 and 27r 
the zeros are splitted in two closely located ones again near the same position in the real 
direction. For lm.9 = 7n/3 and 8n/3 there are triplets of close zeros, etc. In fig.|13|the same 
is examplified near the next zero of T(9) at 9 = 2.2527 .... It is seen already that the scale 
of splitting becomes very small with Re 9 growing and such zero multiplets look like multiple 
zeros if the numerical resolution is not enough. 

3. General real b 2 . In general the periods of T(9) and T(8) are incommensurable. As 
it was mentioned in sect|| it leads in particular to quite complicated Re 6* — > oo asymptotics 
at sufficiently large lm8. While the analytic structure of T{9) and T{9) remains essentially 
as described above (in particular, I verified for many examples that all zeroes of T and T 
are on the real axis), the structure of zeros in X{9) becomes, as lm.9 comes essentially large, 
rather chaotic. I hope to turn again to this point in future studies. Let me mention only the 
following observation conserning the small R (or large negative x) picture. If —x ^> 1, in the 
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Figure 13: The same as in fig. 12 near the next multi-zero at Re# »s 2.2527. . . . 



central region x < 9 < —x function X(9) matches extremely well the following expression 

cos(2QP6) 



X{9) 



(7.3) 



[sinh(27r6P) sinh(27rP/&)] 1/2 

where P is an P-dependent parameter. Roughly it can be estimated from the requirement 
that X(9) ~ at 9 = ±x, i.e., P = ir/(—4Qx) + 0(x~ 2 ). Notice that substituting of this 
approximation to the expression of the effective central charge 

1 - 24P 2 



Ceflf 



(7.4) 



(see ref. |16| for the motivations) we arrive just at the leading UV logarithmic correction 
( p. 161 ). It is easy to verify that expression ( [7.3|) satisfies exactly the functional X-system 
( P- 13] ). This means in particular that it remains a valid approximation (at large —x) of X{9) 
in the whole complex strip x < Re 9 < —x. 



In fact along the lines of [1(J and |L7| a far better estimate of P can be found which takes 
into account all logarithmic in R corrections to Q2.16Q . In this framework P is determined 
as the first root of the transcendental equation 



4QPlog(P/27r) +il g(-S L (P)) = -it 
where Sl(P) is the so-called Liouville reflection amplitude (for the arguments see 

T(b 2 ) \ ~ 2iP/b T(l + 2i6P)r(l + 2iP/b) 



Sl(P) 



7T/i 



r(i-6 2 ) 



T(l - 2i6P)r(l - 2iP/b) 



(7.5) 



(7.6) 
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Figure 14: Function X(9) at x — —6 and different values of b 2 (solid lines) compared with 
the approximation (7.3). Corresponding values of P are indicated. 



For expample, in fig.[14| the shape of X{9) is compared with the approximation ( |7.3| ) for 
x = —6 and two values of the parameter b 2 = 0.8086 . . . and b 2 = 0.1201 . . . According to 



( \l.6\ ) they correspond to P = 0.05332 . . . and P = 0.04052 . . . respectively. 

In view of (|7.3|) the plateau heights of T(9) and T{9) in the central region x < 6 < —x 
and x — > — oo can be estimated as 

T = 2 cosh(27r6P) (7.7) 
f = 2 cosh(27rP/6) 

with the same P determined by ( |7.6| ). 

4. Complex (staircase) values of b 2 . The staircase version of TBA equation (|2.3j ) 
(with the kernel ( |6.3| )) is solved numericaly in the same way as in the ShG case. The structure 
of the solution has some interesting differences from the case of real b. The peculiarities are 
more manifested if parameter 6q in ( |6.3|) is taken sufficiently large and the deep UV region 
— x ^> 1 is considered.. At 9q 3> 1 parameter b is close to imaginary unity and 

is small. 

In fig.[l5| two pictures of X{6) on the real axis are plotted for x = —46 and x = —91, 
both at 8q = 20. The characteristic staircase behavior (in fact very similar to that of 
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Y{9) observed in ref. [ 10 ) is very apparent. It is seen also that expression (|7.3| ) with the 



parameter P determined as the solution to the staircase version of ( |7.5|) (in the present 
case P = 0.06680 . . . and P = 0.04176 . . . respectively) still follows in the central region 
the average behavior of the solution. The deviations (or in other words the corrections to 
the approximation ( |7.3| )) are now oscillating and much bigger then in the ShG case. At 
sufficiently large — x (many amounts of 9q) the ascending part (starting from 9 = x — 9q) of 
the staircase consists of a succession of almost flat steps of constant width 9q, the heights 
being very well fitted by the expression 

_ sm(2QPn/9 ) 
An " sm(2QP/9 ) 1 J 

with n = 0,1,2,.... It should be noted that for the reasons to be explained just below, in 
the case of complex b expression ( |7.3[ ) does not give an approximation in the whole strip 
x < Re 9 < —x. Although (|7.3| ) is still an exact solution to the staircase X-system (|6.6|) , its 
validity is restricted to a certain parallelogram in the complex (9-plane. 



To see this in fig .[16] I present the location of zeros of the function X{9) (at 9 = 20 
and x = —46) in the upper half-plane. The picture is obviously reflected to the lower half- 
plane by the symmetry of X(9). We see several strings of zeros accumulating along the 
lines lm.9 = 7m/2, n = 2,3,4,... which start at the values Re#2 = 47.1500..., Re#3 = 
27.12267. . . w Re 2 - 9 , Re# 4 = 7.11742 . . . w Re# 3 - 9 etc. There are opposite sets of 
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Figure 16: Zeros of X(6) at do = 20 and x = —46. Solid lines stand for very dense sets of 
zeros. 



strings symmetric with respect of the imaginary axis. As n grows they are shifted in steps 
of — #o and #o respectively and meet each other at a certain value of n (n = 5 in the present 
example). After that, continuous lines of dense zeroes extending from — oo to oo are formed. 
In general at large 9 the first such line is already dense enough to produce the effect of 
a "cut" where the behavior of X(9) changes drastically. In our example, across the first 
"cut" at Im9 = 5n/2 the absolute value of X{6) jumps by many (10 5 ) orders of magnitude. 
This is examplified in fig .[II] where X{6) is plotted along the imaginary axis. Notice that 
before lm.9 = 57r/2 function X{9) is almost constant. As 6q grows this effect becomes more 
and more dramatic and in the limit 8q —>■ oo the lines of zeros become real cuts. In the 
forthcoming paper [IS] I'm going to comment more on this effect which plays a crucial role 
in the analytic connection between the staircase behavior at finite 8q and the "sin-Gordon" 
solutions corresponding to purely imaginary b = i/3. 

In fig.[18| location of zeros of the periodic function T{9) is shown at the "moderate" value 
#o = l an d different x. The picture is presented here just to demonstrate two observations: 
(a) At staircase values of b zeros of T{6) are not exactly on the real axis but displaced slightly 
to the complex plane, the displacement becoming negligible very quickly with the number 
of zero, (b) At x — > — oo the picture of zeros is "frozen" in the sense that at x sufficiently 
large the pattern of zeros is simply shifted by the amount Ax in the 9 plane as one changes 
x — > x — Ax. As it is usual in the TBA practice, it is convenient to study these frosen 
limiting patterns substituting the original "massive" TBA equation (|2.3|) by the massless 
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Figure 17: X{6) along the imaginary axis (6* = 20, x = —46). Plot (b) is in the deciloga- 
rithmic scale. 
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version of it. This is one of the motivations for the subsequent study (to be published (18 

8. Concluding remarks 

• In the present study I didn't touch at all the important question about the R (or 



x) dependence of the effective central charge (|2.2|) determined through the finite-size 
ground state energy ( |2.1| ). The UV behavior at x — > — oo is especially interesting 
since the analytic structure of c e g(R) is quite unusual (like in ( |2.16|) ). The Liouville 



quantization condition ( [7.5[ ) together with ( |7.4p proves to be a very good approximation 



to the UV effective central charge behavior. However, while it takes into account all the 
"soft" (logarithmic in R) contrubutions to the asymptotic, there are certainly power- 
like corrections in R. The most important of them (at least at real b) is the contribution 
of the ground state energy ( |1.9|) . Approximation (|7.4j ) is essentially impruved if it is 



taken into account 

c eff = 1 - 24P 2 + 3(mP) 2 /(4tt sin rrp) (8.1) 

Usually after the ground state energy contibution is subtracted the reminder is a series 
in the in the perturbative powers of R like in (|2.15|) . In our present case with Lagrangian 
( |1 . 1| ) naively one could expect a series like 

oo 2n 

c cff - 3(mP) 2 /(4vr sin np) = 1 - 24P 2 + C «( P ) (/^ 2+2&2 ) " (8-2) 

n=l 

where the powers of R are predicted on the dimensional arguments and the coefficients 
c n (P) are something like the Coulomb gaz perturbative integrals corresponding to 
expansion in fi around the vacuum of momentum P (therefore they keep some smooth 
R dependence). Although the leading correction n — 1 of order R 4+4b is likely in 
agreement with the numerical data (at least at sufficiently small b 2 ), the whole structure 
( |3.2|) certainly contradicts duality. Correct expansion must contain also the "dual" 
corrections with powers 

R 4+4/b\ of course many people immediately propose to add 
the dual interaction to the Lagrangian and consider an action like 



.4 



trivial 



^-{d a <p) 2 + 2//cosh(260) + 2/icosh(26" 1 ( 

47T 



d 2 x (8.3) 



with jl taken from ( |1.8|) or sometimes introduced as an independent coupling. Simul- 
taneous expansion in both couplings would supposedly lead to a self-dual series 

oo 

c cfr - 3(mP) 2 /(47r sin Tip) = 1 - 24P 2 + c mtn {P)fi 2m fi 2n R iQ{mb+n/b) (8.4) 

m,n=l 

Coefficients c m>n (P) are computed as the mixed Coulomb gaz integrals which incude 
both kinds of charges produced by the expansion of action ( |3.3| ). Although the general 
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structure of ( |S.4| ) looks very likely, to my sense this scenario (which I call trivial) 
with naive addition of dual interactions (as in ( |3.3|) ) is not exactly in the spirit of 
duality. However at present it does not contradict any data and in fact should be 
verified. A check of this trivial scenario (as well as any other one) requires very subtle 
measurements of the subleading power-like corrections to the effective central charge 
as well as tedious calculations of mixed perturbative integrals. I understand that this 
is a quantitative work which hardly can be replaced by general speculations. 

Of course the definition of P as the solution to the quantization condition (|7.5|) is not 
completely unambiguous. The power-like in R corrections (which are exponentially 
small in 1/P) can be arbitrarily redistributed between the expression for the observable 
effective central charge ( |3.4|) and the formulation of the quantization condition. In other 
words one can add exponentially small in 1/P corrections to ( |7.6| ) and consider this as 
a new definition of P. The problem is that at present the parameter P is not precisely 
observable, i.e., it cannot be directly measured in the TBA calculations (apart from the 
abovementioned definition through the observable finite-size effective central charge). 

At this point we arrive at the most intriguing question touched only slightly in the 
present study. This is about the possibility to construct the solutions Q±{9) and 
Q±{9) of eqs. (|4.11|) with the properties Q4.12D and fl4.13|) such that the solution to 
Hf ) can be built as the combination ( |4.14j ). Was this be possible we'd have another 



unambiguous definition of the parameter P as the Floquet index in ( |4.12|) . However 



there are serious doubts that such solutions exist in any sense, at least for rational 
values of b 2 . To clarify this point, in the next publication [|18| I'll consider the massless 



version of the ShG TBA equation where the parameter P in introduced from the very 
beginning instead of the scale parameter R. In this context the solutions (|4 . 1 2|) can be 
found at least as formal series for irrational values of b 2 in the way that the construction 
( |4.14| ) can be given an exact sense. 



In connection with the periodic structures encoded in the functions T(8) and T(8), it 
seems quite interesting to understand better the R dependence of the coefficients T n 
and T n in the expansions (|4.9|) or (|6]9|) . I checked numerically the R — > asymptotic of 
To and To. The leading asymptotics of To and To (remember that for definiteness in this 
study I always take b < 1, in particular the data discussed in this item were calculated 
at b = 0.3466) are very well fitted by the expressions ( |7.7| ) with P the solution to (|7.5|). 
The correction to ( |7.7| ) for T can be set in the form 

T = 2 cosh(27r6P) + (mP/27r) 4feQ T (1) (P) (8.5) 

where T (P) is a smoothly varying function of P with certain UV limit T (0). As 
for T the correction is better fitted as 

T - 2 cosh(27rP/6) ~ A(P) (mR/2n) a (8.6) 
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with a again numerically close to AbQ. The common TBA experience would expect 
from the periodic structure of T{6) a much smaller correction ~ {mR/2n) A Ql h . A prob- 
able explanation is that there are some power-like corrections to the Liouville quanti- 
zation condition (|7.5| ) so that the "correct" value of P C orrect is off from P (calculated 
from ( |7.5| )) by an amount of order (mR / 2%) 4:b Q 

P = Correct + Pl.O {mR/2nf bQ + ... (8.7) 

With this P C orrect an asymptotic 

f = 2 cosh(27rP corrcct /6) + (mP/27r) 4Q / b T (1) (P corrcct ) + . . . (8.8) 

must hold. If the power AQjb ^> AbQ (like in the present experiment with AQ/b = 
37.31 . . . 3> AbQ = 4.480 . . . ) we can even try to relate the coefficient A in (|8.6|) to the 
leading power correction to ( |7.5|) . As the corrections in (|8.8|) are much more suppressed 
at R — > this seems reasonable. 

It is easy to verify that the leading asymptotic of the higher coefficients T n and T n in 
the expansion ( |4.9| ) are of the form 

T n = (mi?/27r) 2bQ|n| Ti 0) (P) + . . . (8.9) 
T n = (mR/2ir) 2Q ^ /b fi 0) (P) + ... 

with Tn\P) and Tn\P) regular at P = 0. I analysed quantitatively the functions 
Ti°\p) and T^\p). Motivated by the constructions of refs. fT^Ull let us introduce 
slightly rescaled functions 

tl (P) = T} 0) (P) (Z(p)y 2bQ (8.10) 
t^P) =f 1 (0) (P) (Z(p)y 2Q/b 

with Z(p) defined in ( |1.7| ). In table [I] the values of ti(P) and ti(P) are compared with 
the following analytic expressions borrowed from refs. []I1L |1^] where these coefficients 
enter the explicit constructions of the "sin-Gordon" (i.e., related to purely imaginary 
values of b) analogs of T(6>) 

,CFT fpl = 4vr 2 r(l + 26 2 ) 

i ^ > T 2 {b 2 )T{l + b 2 + 2ibP)T{l + b 2 -2ibP) 1 ' J 



^ FT(p); lrT(H 2b--) 



T 2 (b- 2 )T(1 + b- 2 + 2iP/b)Y{\ + b- 2 - 2iP/b) 



The TBA numbers are measured at b = 0.3465545. . . . and different values of x. The 
same for the staircase example 8$ = 10 is presented in tableQ. The convergence 
is noticably slower due to much weeker suppression of the higher power corrections 
(Re(4 + 46 2 ) = 0.1926. . . in this case). The numbers quoted make it clear that an 
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X 


P 


j-TBA 

L l 


fCFT 


rTBA 


ICFT 

L l 





0.3896985 


-0.7228948 


-0.7213854 


-1.452134e-2 


-1.437650e-2 


-2 


0.1152693 


-0.6585087 


-0.6585077 


-8.578078e-3 


-8.578008e-3 


-4 


6.046935e-2 


-0.6542746 


-0.6542746 


-8.272562e-3 


-8.272562e-3 



Table 1: Numerical values of 



X 


P 


J.TBA 
L l 


6 1 




-15 
-30 


0.1771863 

8.604775e-2 

5.695833e-2 


0.0711632-1. 1669319i 

-0.0660054-0.7447625i 

-0.0708064-0.6568688i 


0.0051770-1. 1896763i 

-0.0684458-0.7445886i 

-0.0709316-0.6568002i 



Table 2: Numerical values of 



analytic continuation of the constructions of refs. |14|,|T^] for the ShG or staircase values 



of the parameter is quite relevant. The corresponding "CFT integrable structures" will 
be shown to have a precise relation to the solutions of the massless versions of ShG and 
staircase TBA equations [IS| where the parameter P = P C orrect is fixed by construction. 



I'd like to mention a quite intriguing recent article |2(J where the author arrived at the 
function ( |2.9|) in a rather different context. It appears as the exact wave function of 
the finite-size sinh-Gordon model in special "7-coordinates" which are the sinh-Gordon 
version of the Flaschka-McLaughlin variables (see |2D| for the details). This gives again 
a motivation to continue the study of the analytic structures related to the ShG TBA 
equation. 
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